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We investigate to which extent a time variation of the gravitational constant or other fundamental 
constants affects the best fit of the Hubble diagram of type la supernovae. In particular, we show that 
a slow increase of G in the past, below experimental constraints, can reconcile the SNIa observations 
with an open zero-A universe. 
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Type la supernovae are an exciting tool for cosmology. They are likely to be good standard candles, and can be 
observed down to z ~ 1. Residual systematic effects, like reddening, if-correction, and intrinsic differences in the 
light curve, can be removed by suitable techniques, so that their intrinsic peak magnitude can be estimated with an 
error of less than half a magnitude down to m ~ 24, where different cosmological models begin to be distinguishable. 
Recent work by Perlmutter et al. (1999, hereinafter P99) and Riess et al. (1998) have shown that SNe la are indeed 
i—\ • a powerful test of world models. The most striking results so far, obtained with 42 deep redshift SNe, and 18 low 
redshift SNe (P99; see also Garnavich et al. 1998) are that a zero-A flat universe is ruled out to an extremely high 
confidence level (c.l.), and that the best-fit flat universe is fl m = 0.28 and fl\ = 0.72. An open, zero-A model is also 
ruled out, at 99% c.l.. These results are particularly important because an Q m around 0.3 is in agreement with a 
number of independent observations, especially concerning the mass of clusters. 

The observed SNe at high redshift are roughly half a magnitude fainter than what a flat matter-dominated universe 
predicts. Such a relatively small difference, and the crucial result that depends on it, deserves close examination (see 
^\ . for instance an alternative physical explanation in Goodwin et al. 1999, and a phenomenological one in Drell et al. 

1999). In this paper we investigate whether a slow time dependence of some fundamental constant can affect the 
result. The reason for such an investigation is intuitive: if the peak luminosity L p of SNe la depends on, say, the 
gravitational constant G, then a relation G(z) will show up as a L p {z), thereby modifying the assumption of standard 
candles. In particular, suppose G increases by AG at a certain look-back time t. Suppose then that L p ~ G~ 7 ,7 > 
■ (see below); then the predicted apparent magnitude will be fainter by 

^ ■ . AG(t) 

> : AmG = 2 ^wWg- w 

^ . Let us call such an effect G-correction. If 7 is of order unity, and G at z ~ 0.5, where most deep redshift SNe are, 
was, say, 50% higher than the present value, then the G-correction would be half a magnitude, and could reconcile the 
SN observations with a zero-A cosmology. Naturally, we are thinking that G increases monotonically with look-back 
time, so that the low redshift SNe, as well as the light curve calibration, remain unaffected by the G-correction. 

SNe la are expected to be standard candles because the explosion occurs when a white dwarf accretes matter from 
a companion and reaches the Chandrasekhar mass 

3.1 fhc\ 3/2 

Mc ^ -r TV , (2) 



m' \G, 

where m! is the mass per electron, independently of the progenitor status and of the accretion history (at least in 
the standard model, see e.g. Woosley & Weaver 1986). It is then very likely that some relation between the peak 
luminosity and the Chandrasekhar mass exists, whatever the precise explosion mechanism is (see e.g. Arnett 1982). 
It seems therefore reasonable to test the hypothesis 

L P ~ G-\ (3) 

with 7 > of order unity. 

Actually, any dependence of L p on a fundamental constant, like the nucleon mass, can give an effect similar to the 
G-correction. In fact, a field theory with a coupling of a Brans-Dickc field to gravity can always be transformed into 
a mathematically equivalent theory in which the field couples explicitely to the matter fields, rather than to gravity, 
and therefore the masses are field-dependent. In the following, for definiteness, we focus however on the gravitational 
constant. 
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II. THE G-CORRECTION 



Let us then formulate more precisely the effect we are testing. There are several models which predict a time 
dependence of G, generally based on a Brans-Dicke coupling. Most of them can be simply parameterized as 

G = G (l + tH ) m , (4) 

where Go is the present gravitational constant, t is the look-back time, and Hq ~ 10~ 10 yr~ 1 is the present Hubble 
length for h = 1 (see e.g. Amendola 1999a). The constraints on G can be expressed as follows 



G 

Therefore we can assume 



G| o = |m|i/ 0<a 10- 1 V" 1 - (5) 



\m\ < a. (6) 



The value of a is of the order of 0.1-0.01, depending on the assumptions (see e.g. Guenther et al. 1996). The relation 
between the look-back time r in units of Hq 1 and the redshift is 

T = tH = f [(1 + z') 2 (l + n m z') - z'(2 + z')Q A ] ~ V2 . (7) 

Jo 1 + z 

We can write then L p = L p q(G / Go) -7 and obtain a G-correction 

Am G = 2.5/3 log[l + r(z)], (8) 

where (3 = mj. The maximum value of (3 is of course crucial to our argument. However, it is difficult to determine 
it with some certainty, especially because the SNe la mechanism is still matter of controversy. It is probably safe to 
assume 7 — 1, which implies (3 < 0.1 — 0.01. If the stronger constraint on m holds, then the G-correction is probably 
unable to alter the conclusions of P99 and Garnavich et al. (1998). However, there are at least two possibilities to 
weaken the constraints on (3. First, the dependence of L p on G can be stronger than 7 ~ 1. Second, the G{t) relation 
can be different from the simple power-law (Q) : for instance, G could be an oscillating function of time, as was 
proposed by Morikawa (1991) to explain the periodicity in the pencil-beam galaxy distribution. In this case, tuning 
the oscillating phase and choosing a period of the order of z ~ 0.5 — 0.8, one can have |G/G| very small at present, 
and larger in the past, so as to escape the present constraints. An oscillating G(t) would be naturally produced if the 
Brans-Dicke scalar field oscillates around its potential minimum. For such a behavior, our power law approximation 
is supposed to hold approximately around z ~ 0.5 — 0.8. 

For these reasons, and because in this paper we are interested in testing to what extent a time dependence of the 
fundamental constants affects the SN results, we allow (3 to vary even beyond 0.1. Actually, the SNe la measures can 
be employed also to estimate (3 itself, as an additional parameter to f2 m , as we will show. 

There is actually a small inconsistency with Eq. ([?]) and in the luminosity distance: since we are assuming G to be 
time-dependent, we should use the Brans-Dicke equation, or some variants of a scalar-tensor theory, in place of the 
Friedmann equation. Or, equivalently, if we leave G constant, and vary one of the other fundamental constants, say 
the nucleon mass, we should insert an explicit coupling between matter and the scalar field, and this would change 
the dependence of the matter density with time in the Friedmann equation (see e.g. Wetterich 1995, Amendola 
1999b). However, the tight constraint on G/G implies that the deviation from the Friedmann equation is minimal, 
so that we make only a small error in using it. There is also the possibility that the non-minimally coupled field 
is itself a dynamical cosmological constant, as in quintessence models (Caldwell et al. 1998, Uzan 1999, Chen & 
Kamionkowsky 1999, Baccigalupi et al. 1999, Amendola 1999c). In these models, a new parameter is needed, the 
scalar field effective equation of state w = p/p, and the calculations should take into account explicitely this. Such 
models will be investigated in another paper. 

III. LIKELIHOOD RESULTS 

Once the G-correction is added to the expected apparent magnitude, we obtain the following relation 

m(z) = M + 5logV L {z;fl m ,n A ) + Am G (z-n m ,n A , (3), (9) 
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where we use the notation of P99, in which M is the Hubble-constant-free magnitude zero-point, and T>l the Hubblc- 
constant-frcc luminosity distance. In Fig. 1 we show m(z) for various choices of the parameters, fixing M. to its best 
fit value. Values of of the order of 0.5 or larger seem to explain the data even for A = 0. To be quantitative, we form 
as in P99 a gaussian likelihood function L(Ai, f2 m , CI a, 0) for the four parameters to be estimated. With respect to 
P99 we adopt two simplifications. First, we neglect the dependence on a fifth parameter, called a in P99, the slope of 
the width-luminosity relation, and assume throughout the P99 best value, a = 0.6. As stated in P99, the dependence 
on a has anyway a small effect on the Hubble diagram. The result we get for = 0, extremely close to the original 
P99 results, confirm that this is indeed an acceptable approximation. The second approximation is to neglect the 
correlations among the photometric uncertainties for the high-redshift SNe; as explained in P99, they are small and, 
again, the comparison between our results and the original ones confirms that this approximation is not harmful. 

As in P99, we marginalize the likelihood function over the parameter M (that is, we integrate it over), we assume 
the prior condition fl m > and, finally, we use the catalog of SNe la published in P99, excluding the six SNe rejected 
for the fit C. We are left then with 16 low-redshift and 38 high-redshift SNe. 

We first study how the confidence regions for fl m , J7a depend on 0. In Fig. 2 and 3 we show the confidence regions 
for = 0,0.2 and 0.5. For = we recover the best fit of P99, (ft m ,fi A ) = (0.7, 1.3), while for = 0.5 we obtain 
(0 TO , Qa) = (0.8, 0.9). For a flat universe, the best fit goes from TO = 0.3 when = to Cl m = 0.5 for = 0.5. While 
a zero-A flat universe is still very unlikely even for = 0.5, now a zero-A open universe with Q m < 0.5 is within the 
90% c.L, with fl m = 0.2 the best-fit zero-A universe. Already for — 0.2, an open A = model with enough matter 
to account for the cluster masses, Cl m = 0.3, is within the 90% c.l. 

In Fig. 4 we try instead to estimate the maximum likelihood value for 0, constraining the universe to be flat, 
Qa = 1 — The confidence regions are elongated along the degeneracy line 0d — 2(O m — 0.3). Here we see that to 
reconcile SNe la with flatness and A = one needs > 1, which seems quite difficult to achieve. The best estimate 
is (0, Q m ) = (0.50,0.54) but the maximum is extremely flat along the degeneracy line. Imposing both the constraint 
of a flat universe and Cl m > 0.2, we find that has to be larger than -0.4 (90% c.L). 

IV. CONCLUSIONS 

The inclusion of time dependence in the fundamental constants that determine the SNe la peak luminosity can 
modify in a significant way the best-fit cosmological model. Here we considered the gravitational constant time 
dependence, but other fundamental constants can be used as well. Modeling both the time dependence of G and the 
peak luminosity dependence on G as power laws, the effect on the SNe peak magnitude is proportional to the product 
of the power law exponents. The maximum allowed value of is not very well constrained, since it depends on the 
precise mechanism that powers SNe la. Moreover, an oscillating function G{t) could escape the local constraints. 

We can summarize our results as follows: 

• in a flat universe, > —0.4 if fl m > 0.2. 

• for |/3 1 < 0.2 there is no significant departure from the standard results; 

• for > 0.2 , the P99 confidence regions move by la or more toward smaller A, and a zero-A open universe with 
O m = 0.3 is within the 90% c.l. 

• for > 0.5, a zero-A open universe with fl m = 0.3 is within the 67% c.l. region, 

• for > 1, a flat cosmology without cosmological constant would be reconciled with the SNe la observations. 



We acknowledge useful discussions with Amedeo Tornambe. L. A. also thanks Emanuela Previtera for assistance 
in the data analysis. 
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V. FIGURES 
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FIG. 1. Hubble diagram for the high-redshift SNe of P99. The thick line is the best fit flat universe solution 
(£} m ,fiA) = (0.28,0.72). The thin line corresponds to (fi m ,fiA) = (1,0). The dotted lines are again for (fi m ,fiA) = (1,0) 
but with increasing G-correction, (3 = 0.1, 0.3, 0.5, bottom to top. 
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FIG. 2. One sigma (67%) confidence regions in the plane (fl m , Qa) for (5 = (shaded region) and, decreasing toward smaller 
Q A , for /3 = 0.2 and f3 = 0.5. 




FIG. 3. Confidence regions (67, 90 and 95%) in the plane (O m , Qa) without G-correction (shaded region), and with a f3 = 0.5 
G-correction. The straight lines correspond to a flat universe and a A = universe. 
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FIG. 4. Confidence regions in the plane (/3, Q m ), fixing £7a = 1 — Q. m . The vertical dotted line is the minimum value of Q m 
that is allowed observationally. The horizontal dotted line is the zero G-correction case. The star marks the best fit parameters. 
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